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^ \ Abstract 

^ ' We analyze scenarios in which some flavour of sneutrino is the next-to-lightest supersym- 
■ 

O ' metric particle (NLSP), assuming that the gravitino is the lightest supersymmetric particle 
00 ' 

Q \ (LSP) and provides the cold dark matter. Such scenarios do not arise in the constrained 

supersymmetric extension of the Standard Model (CMSSM) with universal gaugino and 
• 1-^ . 

^ ■ scalar masses input at the GUT scale. However, models with non-universal Higgs masses 
■ (NUHM) do allow scenarios with a sneutrino NLSP, which are quite generic. We illustrate 

how such scenarios may arise, analyze the possible metastable sneutrino lifetime, and ex- 
plore the theoretical, phenomenological, experimental and cosmological constraints on such 
scenarios. We also discuss the collider signatures of such scenarios, how they may be distin- 
guished from neutralino LSP scenarios, and how different flavours of sneutrino NLSP may 
be distinguished. 
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1 Introduction 



In the framework of supersymmetry with conserved R parity, the hghtest supersymmetric 
particle (LSP) can have neither electromagnetic nor strong interactions: otherwise it would 
have bound to conventional matter and been detected in searches for anomalous heavy nuclei 
[1]. Within the minimal supersymmetric extension of the Standard Model (MSSM), weakly- 
interacting candidates for the LSP are the lightest sneutrino u, the lightest neutralino and 
the gravitino G. The (left- 'handed') sneutrino LSP hypothesis is excluded by a combination 
of neutrino counting at LEP and direct dark matter searches [2]. Accordingly, general 
attention is focused on the neutralino (NDM) [3] and gravitino dark matter (GDM) [4, 5] 
possibilities, and in this paper we assume the latter. 

The next question is the possible nature of the next-to-lightest supersymmetric particle 
(NLSP) in a GDM model. Two natural possibilities are the other candidates for the LSP, 
namely the sneutrino and the neutralino, but the NLSP could equally well be charged and 
even coloured. Indeed, the lighter stau slepton is a natural candidate for the NLSP [6] 
within the constrained MSSM (CMSSM) with gravity-mediated supersymmetry breaking, 
in which the soft super symmetry-breaking scalar masses mo, trilinear parameters Aq and 
gaugino masses mi/2 are each assumed to be universal at a GUT input scale [7,8]. The 
lighter stau is also a natural possibility within minimal supergravity (mSUGRA), in which 
the gravitino mass is fixed: rriQ = rriQ, and there is an additional relation between trilinear 
and bilinear soft supersymmetry-breaking parameters [9]. Another possibility for the NLSP 
within a scenario with non-universal Higgs masses (NUHM) [10-14] is the lighter stop [15]. 

In this paper we study the other possibility for the NLSP within the NUHM with gravity- 
mediated supersymmetry breaking, namely the lightest sneutrino, assuming that the grav- 
itino provides the cold dark matter 0. Like other NLSP candidates in gravity-mediated 
scenarios, the sneutrino NLSP within the NUHM is expected to be very long-lived. The 
dominant decays of the other NUHM NLSP candidates produce particles with copious in- 
teractions such as charged particles and photons, which are subject to strong cosmological 
limits [17]. These limits are strong enough to exclude effectively all of the parameter space 
where the lightest neutralino is the NLSP [18] [^. However, the dominant decay mode of the 
sneutrino is z/ — Gu, and the cosmological limits on neutrino injection are much weaker than 
those on the injection of photons and charged particles [19,20]. Therefore the cosmological 
limits on NUHM z> NLSP scenarios are relatively weak, leaving considerable scope for a sneu- 

^The possibility of a sneutrino NLSP has also been studied within gaugino-mediated models of super- 
symmetry breaking [16]. 

^If the gravitino mass is much lighter than the neutralino mass, these limits might still be satisfied. 
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trino NLSP. On the other hand, the sneutrino must appear at the end of the decay chain of 
every MSSM sparticle produced at a colhder, and the particles produced in supersymmetric 
decay cascades provide distinctive experimental signatures for a sneutrino NLSP [21]. In 
particular, the charged leptons produced in association with a z/ NLSP provide tools for 
diagnosing its flavour. 

The structure of this paper is as follows. In Section 2 we discuss sneutrino properties 
within the NUHM, including its mass and lifetime. In Section 3 we discuss the relic abun- 
dance of sneutrinos after freeze-out from a primordial plasma in thermal equilibrium. In 
Section 4 we analyze the NUHM parameter space and identify regions where the NLSP may 
be either the i>e,^ or the v-r- In Section 5 we discuss the cosmological constraints on i) NLSP 
scenarios, and show that they are not severe. In Section 6 we discuss some signatures of 
metastable sneutrinos in different NUHM scenarios, in particular those with different lepton 
flavours accompanying the z/ NLSP. Our conclusions are summarized in Section 7. Calcula- 
tions of the three-body decay v Gwj are described in an Appendix. 

2 Sneutrino NLSP Properties 

We assume that the GUT scale is the effective input scale at which the soft masses are 
specifled, presumably via some gravity-mediated mechanism, and make the NUHM assump- 
tions that the gaugino masses are universal, as are the squark and slepton masses, whereas 
the soft supersymmetry-breaking contributions to the Higgs scalar masses are non-universal. 
We then calculate the physical supersymmetric mass parameters at a low energy scale from 
the running given by the renormalization-group equations (RGEs). We assume that the 
right-handed neutrino supermultiplets, being singlets, get very large Majorana masses, and 
therefore decouple from the low-energy effective theory □. The sneutrino NLSP discussed in 
this paper is essentially the scalar partner of some left-handed neutrino. The flavour of the 
9 NLSP is, however, model-dependent, as we discuss below. 

^However, if the decoupling energy scale is significantly below the GUT scale, the low-energy spectrum 
may be affected [22] through the RGEs. This effect is neglected here. 
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2.1 Sneutrino Mass 

In order to calculate the sneutrino mass, we first look at the RGEs of the slepton sector 
[23-25]: 

dt 87r2^ 2 ^' 

1 

^ = 5^(-49?Mf + 4S). 
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-3^/2m| - glM', + hl{m'~ + + + Al) - 2S), 



^ = ^(-4^7?M,^ + 2/.^(m|^^+mf^ + m? + A?)+45), 



where 



O 9l / 2 2io^2 2 o2| 2, 2\ 

b = — {rrir, — mi + 2{m~ — m~ — 2m~ + m~ + m~ 

+ (m~ -m~ -2m~ + m| +m~)). (2) 

^ QSL L-iL tR bR ' TR>) V I 

The term vanishes and does not contribute in models with universal soft masses such as 
the CMSSM. However, for non-universal models, S can be far from zero and its contribution 
to the RGEs can be significant. When S* = 0, and assuming universal soft breaking masses 
for the L and R sleptons at the GUT scale, the R slepton is lighter than the L slepton at 
the weak scale. However, since S contributes in opposite ways for L and R sleptons, if S is 
large and negative, could be lighter than L^. Furthermore, there are additional D-terms, 



2 2 

m.^ = rrij 

2 



— cos(2/5)m| — sin^ 9w 



m-, 



- m|^+cos(2/5)im|> (3) 

that split the sneutrino and charged-slepton masses. Since cos(2/5) < for tan/3 > 1, the 
sneutrino is lighter than its charged-slepton partner. 

Comparing the first two generations and the third generation, we see that the tau sneu- 
trino could be lighter than the electron and muon sneutrinos because of the Yukawa terms 
in the RGEs. However, this is not always the case, due to the fact that m\ could be negative 
and large, and we display examples later where the z/ NLSP has electron or muon flavour. 
The lighter stau mass is also suppressed by off-diagonal terms in the mass matrix. Thus, 
depending on the model parameters, either the tau sneutrino or the lighter stau might be 
lighter. 



''Although we write the one-loop RGEs for simplicity, our calculations include two-loop contributions. 
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Therefore, a sneutrino could be the NLSP if S is large and negative @. We see from ([2]) 
that 5* is negative when m\ — m\ < o|^, which is not possible in the CMSSM, for which 5 = 
by assumption. We now study how this may occur in the NUHM model, using the freedom 
that the Higgs soft supersymmetry-breaking masses at the GUT scale are not necessarily 
equal to mo, the universal scalar mass for sleptons and squarks. 

The electroweak symmetry breaking conditions may be written in the form: 

m\{Q) = m?(Q) + ml{Q) + 2ii\Q) + A^(g) (4) 

and 

„ m? -m^tan2/5+ (1 -tan2/3) + AW 

= — — ' 5 

tan2/?-l + A(f) 

where A^ and A^j^'^^ are loop corrections [23,24,26] and mi 2 = Tni 2{m^z) B The values of 
the NUHM parameters at Q are related to their values at m^^ through the known radiative 
corrections [23,25,27] ci,C2 and c^: 

m\{Q) = ml + ci, 
ml{Q) = m\ + C2-, 

= /i' + c^- (6) 

Solving for m\ and m^, one has 

mi(l + tan2/3) = m^(Q) tan^ /? - /i^(tan2 /5 + 1 - Ajf^) - (ci + ca + 2c/,) tan^ /? 

1 
2 

and 



-Aa(Q) tan2/5 - -mlz{l - tan^/?) - A« (7) 



ml{l + tan^ P) = m^(Q) - fi\tan^ /? + 1 + AJ^^)) _ ^ ^ 2c^) 

-Aa{Q) + \ml,{l - tan^/?) + A« . (8) 

The correction A^^^^ is positive and generally of 0(0. 1). From here we can see that there are 
two possible ways to get negative S via negative m2 — m^: the first is by using very large /i^, 
and the second is by using very large m\. If m\ is relatively small while /i^ is very large. 



^One could also obtain a light sneutrino within a supersymmetric SU(5) GUT with different soft masses 
for the 10 and 5 multiplets [20]. Another alternative is within a gaugino- mediated supersymmetry breaking 
model, in which the Higgs masses are again different from the other sfermion masses [16]. 

^Assuming that these are the dominant terms in 5, which is the case for the NUHM that we consider 
here. 

^Our convention is such that Hi^2 = Hd^u, and tan/3 = V2/V1. 
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the Higgs masses-squared mjj_^ ^ = 2 + may be negative at the GUT scale, which could 
lead to a vacuum stability problem [28] . 

A weak-scale scalar mass-squared with m^{MGUT) < generally produces a vev of order 
the weak scale that disappears as the RGEs are run down to the weak scale. Such negative 
masses-squared are not dangerous. It may happen, however, that an instability occurs along 
some F— and D— fiat direction. In this case, a negative mass-squared may be large and 



lich the Universe becomes 
This possibility has been 



still present at a renormalization scale, Q ^ v [29]. Models in w 
trapped in such non-Standard Model vacua are clearly excluded!^ 
studied in the NUHM along the Hi - H2 and H2 - L fiat directions of the MSSM [31]. We 
delineate below regions in the parameter plane where these vacua may be problematic. 

2.2 Sneutrino Lifetime 

In the GDM scenario used here, the sneutrino NLSP would eventually decay into the grav- 
itino, and the dominant decay channel is the two-body decay 

u-^G + u, (9) 

with the decay rate 



2b 



487r Ml, \ m? 



'PI' G 




(10) 



where m~ is the gravitino mass and Mp\ is the Planck mass: Mpi = l/y^SvrGjv — 2.4 x 
10^8 GeV. 

We plot in Fig. [T]the sneutrino lifetime, ^ l/r2b, as a function of the gravitino mass 
for rric, = 10, 100, 500, and 1000 GeV respectively. Note that we plot the lifetime only for 
Am = m^; — > 1 GeV. Clearly, a smaller mass gap would yield an even longer lifetime. 
We see that the sneutrino lifetime could be less than 1 second only when mc, is large, or the 
gravitino mass is (much) less than 1 GeV. On the other hand, if the mass gap is small, the 
sneutrino lifetime can be very long, potentially even longer than the age of the Universe, 
which is 0(10^'') s. However, there are cosmological and astrophysical constraints on the 
possibility of a sneutrino with lifetime longer than the age of the Universe at recombination 
that we discuss in more detail later. 



^Whether this occurs or not depends on the specific cosmological history during inflation [30,31]. 
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Figure 1: The sneutrino NLSP lifetime as a function of for rrij^ = 10,100,500 and 
1000 GeV (top to bottom). 

3 Cosmological Sneutrino Density 

Assuming thermal equilibrium in the early Universe, one can calculate the sneutrino relic 
density after decoupling but before its decay into the gravitino. This is done by the usual use 
of the Boltzmann equation and calculation of the sneutrino annihilation and coannihilation 
cross sections. The calculations are identical to those required to calculate the relic sneutrino 
density if it is the LSP. 

The possible sneutrino-pair annihilation two-body final states are the following [12]: 



Initial State 


Final States 




//, W+W-, ZZ, hZ, hA, HZ, HA, hh, hH, HH, AA, AZ, H+H', W+H', H+W- 
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If the soft masses for the sfermions are universal as assumed here, the electron sneutrino 
is always degenerate with the muon sneutrino, and the tau-sneutrino mass might be nearby. 
In the NUHM case that we consider here, there could also be other sparticles that are almost 
degenerate with the sneutrinos, such as the lightest neutralino and chargino, and charged 
sleptons. We list below the corresponding coannihilation processes and their possible final 
states: 

• Coannihilation with other sneutrino flavours: 



Initial State 


Final States 











• Coannihilation with charged sleptons: 



Initial State 


Final States 




//', hE-, HH~, AH-, hW-, HW-, ZW', -fW, W~A 











In the first line above, the coannihilation is between a sneutrino and a charged slepton 
of the same generation, whereas in the second and third lines they are not necessarily 
from the same generation. 

• Coannihilation with the lightest neutralino: 



Initial State 


Final States 




uZ, vh, vH, vA, e-W+J-H+ 



and similarly for u*. 
• Coannihilation with the lighter chargino: 



Initial State 


Final States 




^7, ez, eh, iH, iA, vW-, vH- 


X'^i^ 


uW+,iyH+ 
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and similarly for u*. 

In our scenario, the sneutrino NLSP eventually decays into the gravitino before the 
current epoch. Consequently, the gravitino relic density is related to the sneutrino density 
before its decay by 

Q h' = !!^Q.h^ + Qlh', (11) 

m^, ^ 

where is the contribution to the gravitino density from thermal production after reheat- 
ing, which is sensitive to the unknown reheating temperature Tr. We do not discuss this 
contribution here. The only constraint we impose is that the contribution to the gravitino 
relic density arising from sneutrino decay does not exceed the value suggested by WMAP [32] 
and other observations: 

(]j3M/i^ = 0.1099 ±0.0062. (12) 

Hence, we require that the first term on the right-hand side of ffTTl) should not exceed 
~ 0.1223 (the 2-a upper limit). Because of the scaling by the mass ratio m-^/my, even 
a large sneutrino density after decoupling could still be compatible with the dark matter 
constraint if mg <^ niy. 

In this case, we must check whether gravitinos are non-relativistic at the time structure 
formation begins, roughly at ~ 5 x 10^^ s, in which case they act in the same way as 
conventional cold dark matter. If m3/2 -C m/>, -E'3/2/'^3/2 — fTiv/'^fTi^/2 and i?3/2 scales 
subsequently as (rj>/t)^/^, where is the sneutrino lifetime. We can use Fig. [T]to estimate 
whether or not gravitinos will behave as cold dark matter. For example, at the left-most 
point of the lowest curve, gravitinos are produced relativistically with -E3/2/'n^3/2 ~ 1000, but 
the decay occurs so early that they become non-relativistic well before structure formation 
begins. The same is also true for the two middle curves in Fig. [TJ Only for the left part of 
the topmost curve (when m/^ = 10 GeV) is there a potential problem. However, in this case 
the sneutrino relic density is generally already small in the model considered here. 

4 NUHM Parameter Space 

In the CMSSM, the values of itla and yU are determined by the electroweak vacuum conditions 
for any given input values of mi/2, ""^o? and tan (3. However, the constraints on tua and 
are relaxed in the NUHM, with the values of these parameters being related to the degrees 
of non-universality assumed for the Higgs soft masses m\ and m\. A general discussion of 
the parameter space of the NUHM was given in [11-14], which we use here as a starting 
point for our discussion. Regions of the NUHM parameter space in which a sneutrino is the 
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lightest spartner of any Standard Model particle were identified in [12], see for example the 
dark blue shaded regions in Figs. 2, 3, 4, 6, 8 and 9 of that paper. There, it was assumed 
that the LSP is the lightest neutralino Xi with the gravitino assumed to be heavy, so that 
these light-sneutrino regions were disallowed. However, in this paper we assume that the 
gravitino is the LSP, so the viability of these light-sneutrino regions must be re-evaluated. 
We focus our discussion here on (/i,m^) planes of the types shown in Figs. 4 and 6 of [12], 
where the light-sneutrino region moves to lower as increases. 

In general, keeping some sneutrino species light favours small values of mi/2 and ttiq. 
However, there is an important lower limit on mi/2, in particular, due to the LEP lower 
bound on the mass of the lightest MSSM Higgs boson. A sneutrino NLSP region may be 
found by choosing a moderate value of mi/2 = 500 GeV while keeping mo relatively small, 
e.g., mo = 100 GeV. The resulting masses of the sparticles are shown in Fig. [2] assuming 
tan/3 = 10 and Aq = 0, for various values of m^ = 200, 1000, 1500 and 2000 GeV, in panels 
(a), (b), (c), and (d), respectively!^. 

As noted earlier, in regions of the NUHM parameter plane (particularly when fi and m^ 
are large), the masses-squared of the Higgs and left-sleptons have a tendency to run down 
to negative values at the GUT scale. This allows for the possibility that large scale vevs be 
excited along the Hi — H2 or H2 — L fiat directions. One expects that these fiat directions 
are lifted by some effective operator at or above the GUT scale. The vev along the flat 
direction is sensitive to the fundamental scale associated with this operator and clearly grows 
as that scale is increased above the GUT scale. The reliability of this high-scale vacuum 
depends also on the one-loop corrections to the scalar potential. This sensitivity can be 
characterized by the ratio of the tree-level vev to the renormalization scale, Qo, at which the 
vev disappears (i.e., the masses- squared go positive). Here we adopt the most conservative 
set of assumptions, namely that the flat directions are lifted at the GUT scale and that the 
vev must be of order Qq {e = 1 in the notation of [31]). This preserves the largest volume of 
the NUHM parameter space. Of course the GUT constraint itself is cosmology-dependent, 
and may not be important if the Universe starts out in the weak scale vacuum after inflation. 
For more details on this constraint see [31]. 

In the case of small m^ = 200 GeV, shown in panel (a) of Fig. [2|, we see that the GUT 
stability constraint (represented by a couple of vertical black dotted lines) allows only small 

^1.1 TeV, far from the sneutrino NLSP region which appears when |/i| ^ 2800 GeV. We 
also see that the 6 57 constraint (solid green line) allows only fi ^ 1600 GeV and the Higgs 
mass constraint (dashed red) allows only /i ^ —1000 GeV for this small value of m^. Note 

^We use mt — 172.6 GeV for our analysis [33]. 
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Figure 2: Sparticle masses as functions of n for tan (3 = 10, mi/2 = 500 GeV, ttlq = 100 GeV, 
Ao = 0, mt = 172.6 GeV, mfe(mfe)^ = 4.25 GeV, and tua = (a) 200 GeV, (h) 1000 GeV, 
(c) 1500 GeV and (d) 2000 GeV, respectively. In panels (c) and (d), the sparticle lines are 
truncated at larger |/i| where some sneutrino becomes tachyonic. Gonstraints are represented 
by vertical lines: black dotted for the GUT constraint (larger is excluded); red dot-dashed 
shows the Higgs mass contour at rrih = 114.4 GeV, while the constraint using the LEP 
likelihood function convolved with theoretical uncertainties in the Higgs mass (computed here 
using FeynHiggs [34]) is shown by the red dashed line; the {g — 2)^ constraint (described 
in the text) is shown by the light blue long dashed lines; and solid green for the b ^ s'j 
constraint (smaller fi is excluded). 
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that a Higgs mass of 114.4 GeV is found around fi ^ —300 GeV (dot-dashed vertical hne) 
for this value of ttla, but current theoretical and experimental uncertainties only provide for 
the weaker bound shown by the dashed line. The anomalous magnetic moment of the muon, 
{g — 2)^, is reconciled with experiment at the 95% CL for fi ^ 1300 GeV, as shown by the 
light blue long-dashed line Therefore, we do not have an allowed region in panel (a), 
assuming the GUT constraint holds. For larger rriA, these constraints become more relaxed, 
and an allowed sneutrino NLSP region emerges. 

For rriA = 1000 GeV, shown in panel (b) of Fig. [2l the GUT stability constraint allows 

^ 1.5 TeV, the b —>■ s'j constraint allows ^ 0, and the Higgs constraint is essentially 
unimportant, as a Higgs mass greater than 114.4 GeV occurs at /i ^ —700 GeV. The — 
constraint is satisfied in two regions: < /i < 250 GeV and fi ^ 1900 GeV. While the GUT 
constraint is relaxed, the sneutrino LSP region still requires 2.4 TeV ^ |/i|, where we have 
degenerate i>e,ij. NLSPs. We note that around fi = 2.5 TeV several other sparticles are only 
slightly heavier than the i>e,^, including the lightest neutralino the and fiL, the lighter 
f and the Ur- Thus, all these sparticles must be included in coannihilation calculations of 
the Ue,^ abundance. For larger /i, only and /i^ masses stay close to the NLSP mass, while 
the others get larger mass gaps. 

When rriA is increased to 1500 GeV, as shown in panel (c) of Fig. [21 the GUT constraint 
remains at ^ 1.5 TeV, very close to the region when the sneutrino is the NLSP, which 
extends from /i ~ 1.6 TeV to ^ 2.5 TeV. The lightest sparticles are again the z/g^^, with 
the X;^^L,/Wi,^i and Ur again slightly heavier. In this case, a theoretical upper limit on |/i| 
arises when the i'e,fi become tachyonic In this case, the Higgs constraint requires only 
that fi > -1300 GeV, and {g - 2)^ is satisfied when < < 400 GeV or fi > 1700 GeV. 
The constraint from 6 ^ 57 is unimportant at this value of m^. 

Finally, in panel (d) of Fig. [2] we display the sparticle masses for rriA = 2000 GeV. In this 
case, the allowed sneutrino NLSP is the Ur, for 200 GeV ^ \^\ ^ 1.1 TeV. The x becomes 
the NLSP for smaller in the Higgsino region, and upper limits are provided by the LEP 
lower limits discussed above. For this value of m^, the differences in mass between the Ur 
and the heavier sparticles are relatively large. Neither the Higgs mass nor h provide a 

constraint, while [g — 2)^ requires < yU ^ 600 GeV or /i ^ 1000 GeV with the later has mj> 
less than the LEP limit. 

In these plots, we find that sneutrino NLSP has relic density of order O(10~^) which 



""^^We assume that the deviation of {g — 2)^/2 from the standard model is between 10.7 to 44.3 x 10^^", 
the 2 (T range according to [35]. 

^^Here and in panel (d), we truncate all the other sparticle lines at this boundary of the tachyonic region. 
Slightly more stringent upper limits on come from the lower limits on mo provided by LEP [36,37]. 
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is well below the WMAP limit. This means that most of the gravitino dark matter must 
be produced by some other sources, e.g., by reheating. We take coannihilation effects into 
account for the relic density calculation. However, in contrast with the neutralino LSP 
case, coannihilations do not always reduce the final relic density. Sneutrino coannihilation 
with the lightest neutralino would indeed generally increase the relic density, while that 
with charged sleptons might reduce it. In the former (latter) case, the effective sneutrino 
cross section is averaged with the weaker (stronger) annihilation cross section of neutralinos 
(charged sleptons). The relatively small relic density of the sneutrino compared to that of 
the neutralino can be attributed generically to the fact that the sneutrino is a scalar particle, 
rather than a Majorana fermion. In Fig. [2](d), for example, since the mass gaps with other 
sparticles are relatively large, the coannihilation effects are not maximal, but the relic density 
is still small. 

To get a more comprehensive view of the NUHM parameter space, we display in Fig. [3] 
some contour plots in selected (/i, m^) planes. In panel (a), we choose tan/3 = 10, mi/2 = 
500 GeV, mo = 100 GeV and Aq = 0. This panel therefore includes and extends the specific 
examples shown in Fig. [21 We plot the regions where the lighter stau, the right selectron, 
some sneutrino or the lightest neutralino is the NLSP. At large the electron-sneutrino 
is the NSLP and is shown by the regions shaded dark blue. At large and smaller 
the NLSP becomes the tau-sneutrino (shaded light blue). Above these regions, the white 
area corresponds to an unphysical region where one or more of the sparticles has a negative 
mass-squared at the weak scale. Below these regions, in most of the area, it is the lightest 
neutralino which is the NLSP. At lower and relatively small \fi\, we see regions where 
the lighter stau (shaded brick red), or the right selectron (shaded orange) is the NLSP. 
At even smaller \fi\ the NLSP is a higgsino-like neutralino. The narrow turquoise shaded 
region is that in which I^nlsp^^ is within two a of the WMAP value. (Recall that, with the 
gravitino as the LSP, this is not the dark matter relic density.) Regions surrounded by the 
strip (which have higher ^nlsp^^) niight also be permitted if mg -C mNLSP, and regions not 
surrounded by the strip (which have lower ^nlsp^^) would certainly be permitted by the 
dark matter constraint. As one can see these regions track very closely the degeneracy lines 
between the neutralino and one of the four sparticles where the relic density is controlled by 
coannihilations or the funnel region where 2m^ ~ niA (The contour 2m^ = is shown by 
the thin blue line.). 

The Higgs mass contour of 114.4 GeV (red dot-dashed line) excludes ^ ^ —700 GeV. 
However, if one uses a likelihood analysis for the Higgs mass and allows for theoretical 
uncertainties in the calculation of rrih, this constraint is relaxed and most of the displayed 
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Figure 3: Some (/i, myi) planes in the NUHM for (a) tan j3 = 10, mi/2 = 500 GeV, itlq = 
100 GeV, and Ao = 0; (h) tan/3 = 10, mi/2 = 500 GeV, uiq = 100 GeV, and Aq = 1000; (c) 
tan/3 = 10, mi/2 = 500 GeV, mo = 300 GeV, and Aq = 0; (d) tan/3 = 40, mi/2 = 500 GeV, 
mo = 100 GeV, and Aq = 0. In each case, we used rrit = 172.6 GeV and mb{mb)^^ = 
4.25 GeV. Gontours and shading are described in the text. 
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region is allowed, as shown by red dashed line. The b ^ S'j constraint (shaded green) excludes 
small rriA and prefers positive fi. However, this constraint is unimportant in this panel at 
large rriA. The GUT constraint is represented by a black dotted line: it excludes large 
and essentially all of the electron-sneutrino NLSP region for this set of parameters. Finally, 
the light pink shaded region bordered by a black solid line represents the region favored by 
{g — 2)^. The vertical dashed black lines correspond to the chargino mass contours of 104 
GeV and exclude very small values of 

In panel (b), we display a case with Aq = 1000 GeV, with the other parameters chosen 
to be the same as in (a). We see that the stau NLSP region becomes bigger and there is no 
longer a right-selectron NLSP region. In the white region interior to the stau NLSP region, 
the stau has gone tachyonic at the weak scale, which is problematic. We also see that the 
b ^ s'y constraint becomes stronger, especially for negative fi. However, the qualitative 
features of the sneutrino NLSP band at large and/or rriA are similar. This feature is also 
retained in panel (c), where the larger value mo = 300 GeV is chosen, and also in panel 
(d), in which a larger value tan/? = 40 is chosen. In panel (c) (which now extends to higher 
values of itla), there is no longer a charged NLSP. In panel (d), we see again a region with 
a right selectron NLSP. The GUT constraint in this case constrains only the lower right 
corner of the plane shown, and allows part of the electron-sneutrino region. The fi^fi^ 
constraint (orange dashed line) excludes a region with small rriA [38]. 

We conclude that the possibility of a sneutrino NLSP is quite generic in the NUHM, and 
certainly not much less plausible than the lighter stau. This is in contrast to the CMSSM, 
where a stau NLSP is a generic feature at large mi/2 and small ttiq, but there is no possibility 
of a sneutrino NLSP. 

5 Cosmological Constraints on a Sneutrino NLSP 

The cosmological impact of a long-lived sneutrino depends on its lifetime 
decays during or after Big-Bang Nucleosynthesis (BBN), it could alter predictions for the 
light-element abundances. If the sneutrino decays around or after the time of recombination, 
it could distort the blackbody spectrum of the CMB. If the sneutrino decays at a very late 
time, its effect might be seen on the diffuse neutrino and photon spectra. The production 
of relativistic neutrinos by sneutrino decays could also change the equation of state and 
therefore the evolution history of the Universe [39]. 

Although the sneutrino is neutral, and its dominant two-body decay channel produces 

""^^See [19] and references therein. 



. If the sneutrino 
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only a neutrino and the gravitino, which are also neutral, there could still be a significant 
effect on BBN if the sneutrino decays during or round after the time of BBN [19,20,40,41]. 
If the mass gap is sufficiently large, the decay of the sneutrino produces high-energy, non- 
thermal neutrinos. Through scattering processes with the background particles, such as 
z/j + z^j,BG ~^ (e^,/i^,T^), Vi + z/j^BG ^ 7r+ + TT^ and z/j + ef^Q ^ 7r° + vr^, the energetic 
neutrinos transfer some parts of their energies to charged particles. The final-state particles 
may then photodissociate or hadrodissociate the elements already produced by standard 
BBN processes. In the case of the charged pion, it can alter the neutron-to-proton ratio 
if it occurs at the beginning of BBN. There can also be energy transfer through elastic 
scattering with electrons and positrons: z^j + —>■ z/j + e^, and the high-energy might 
then initiate electromagnetic showers. However, at the epoch of interest (when their energies 
are 0{1) MeV) the electron and positron number densities are already low. Therefore these 
processes can be neglected. 

Subdominant three- and four-body sneutrino decay channels can also be important, even 
though their branching ratios are relatively small. This is because these decays produce 
charged and/or strongly-interacting particles directly. These effects had been studied in [19], 
where it was found that the effects of the three- and four-body decays are negligible if their 
collective branching ratio is less than about 10~^. To estimate this branching ratio, we 
calculate the decay rate for the following process 

V~^G + iy + -f, (13) 

which occurs through neutralinos exchange, with the photon produced via the photino con- 
tent of the neutralinos. This provides an estimate of the total multi-body decay rate that 
should be accurate to within an order of magnitude. The detailed calculation can be found 
in the Appendix. 

To gauge the possible impact of the BBN constraint, we examine the tau sneutrino NLSP 
region in Fig. [2](d), where the sneutrino mass varies from the LEP lower limit of about 
40 GeV up to about 100 GeV. We first consider the case mi/2 = 500 GeV, rriQ = 100 GeV, 
tan/? = 10, Aq = and = 2000 GeV shown in the top two panels of Fig. HI Panel 
(a) shows the three-body-decay branching ratio for various values of mg = 1,10,30 GeV, 
corresponding to u lifetimes ^ 10^ s. We see that the branching ratio is always very small, 
falling below 10~^ throughout the range of parameter space considered. This is consistent 
with the results of [20] which also finds a small hadronic fraction when the sneutrino mass 
is ^ 100 GeV. Thus, according to the analysis of [19], the three-body V decay is too small 
to affect significantly the successful results of BBN. 
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tan |3 = 40 , = 500 , mo = 100 , Aq = tan p = 40 , mi/2 = 500 , = 100 , Aq = 




uiv (GeV) my (GeV) 



Figure 4: Exploration of the BBN constraints on sample tau-sneutrino NLSP points with 
/i > 0; (top) mi/2 = 500 GeV, mo = 100 GeV, tan (3 = 10, Aq = and tua = 2000 GeV 
[cf. Fig. W(d)] and (bottom) mi/2 = 500 GeV, mo = 100 GeV, tan/5 = 40, Ao = and 
mA = 1300 GeV [cf, Fig.\3(d)]. Panels (a, c) display the three-body- decay branching ratios, 
and panels (b, d) the i^(T,e) '^'e.lic density. 
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Although the hadronic branching ratio is expected to be small, BBN nevertheless sets a 
limit on the density of sneutrinos at the time of decay. For a 100 GeV sneutrino, all lifetimes 
are safe so long as the quantity 

y^M^ = n^h^ X (3.65 X 10'^ GeV) (14) 

is less than 0(10"^^) GeV for Bh = 10'^ and less than 0(10"®) for Bh = 10'^ where is 
the ratio of the number density of sneutrinos to entropy, nj^/s. In panel (b) of Fig. IH we 
display the sneutrino relic density following freeze-out but prior to decay as Yf^Mj^. We see 
that YuMc is always below about 10^^^ GeV, with a large dip at ~ 45 GeV due to the 
Z resonance in sneutrino-pair annihilation (with a smaller dip at m,> ~ 60 GeV due to the 
h resonance). Thus, the sneutrino density is far below the range where BBN constraints 
become important for the range of the three-body branching ratio shown in panel (a) of 
Fig.H 

The lower panels of Fig. H] display the three-body branching ratio and YaMj;, for another 
case: mi/2 = 500 GeV, mo = 100 GeV, tan/5 = 40, Aq = and rriA = 1300 GeV. We see in 
panel (c) that the three-body branching ratio is smaller than 10~^ for < 110 GeV, and 
always < 3 x 10~^. Panel (d) shows that F,>M,> is, again, always below about 10^^^ GeV. In 
this case, the prominent dip due to the light h resonance in sneutrino-pair annihilation, and 
the direct-channel Z resonance is less important. This difference from the previous case is 
due to the larger value of tan /3. 

These examples are indicative that the sneutrino LSP regions in the NUHM parameter 
space are generally safe from BBN constraints. We next check other possible constraints on 
decaying sneutrinos. 

When the high-energy decay neutrinos are thermalized, their energy is transferred and 
converted to radiation. If the sneutrinos decay after about z = 10^ (corresponding to a 
lifetime of 1.8 x 10^ s), then the photons produced might not have a chance to thermalize, and 
could show up as distortion of the CMB black-body radiation spectrum. These constraints 
were also considered in [19] where it was found that for lifetimes between 10^ and lO^^s, 
the upper limit on Yt^Mp is roughly 10"^ (10"^) GeV for Bh = 10^3(10"^). Thus in the 
parameter space we are interested in, this too is never a serious constraint. 

When the sneutrino and gravitino mass are nearly degenerate, the sneutrino lifetime 
might be very long (^ 10^^ s as seen in Fig. [1]). If the sneutrinos have decayed after the time 
of recombination, they will produce a diffuse neutrino and photon background. In principle, 
there is then a very strong constraint from water Cerenkov detectors placing an upper limit 
on YyMc, of order 10^^^ — 10^^^ GeV [40]. However, these detectors lose sensitivity for 
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neutrino energies below about 10 GeV [42]. Thus we again do not expect a severe constraint 
placed on the parameter space of interest. 

6 Signatures of Metastable Sneutrinos at Colliders 

We have seen in the previous sections that a sneutrino NLSP is a generic possibility in the 
NUHM. It would be metastable, so that its decays would not be seen at colliders, but the late 
decays of relic sneutrinos are not excluded by the available cosmological constraints. The 
mass of such a sneutrino NLSP might be as low as the LEP lower limit [36,37]. As a non- 
decaying neutral particle, the sneutrino would have a missing-energy signature at colliders. 
Distinguishing the sneutrino from other possible origins of such events would require a search 
for the heavier states that decay into the sneutrino inside the detector. 

Covi and Kraml [21] have studied several scenarios with a sneutrino NLSP assuming 
the following mass hierarchies, (a) rrif > > rrii)^: in this case, neutralino decays into a 
neutrino and sneutrino are invisible, and the signatures of such decay chains resemble those in 
a conventional neutralino LSP scenario. However, if > nii^^+mr, the neutralino can decay 
into X ~^ 'T^rf f 1 where the //' pair is soft if the mass gap is small, (b) > nif > irty^: in 
this case, besides the invisible v + i) decay mode, the neutralino can undergo cascade decays 
that might be detectable, such as r + f. (c) > rngj^ > nif > rriy^: in this case, there 
are additional decay channels with neutralino decay into an electron and a selectron which 
subsequently decays into leptons and a tau-sneutrino. Heavier sparticles might decay via 
the lightest neutralino, but they might also decay directly into the e.g., via X2 ^ ^ + 
or via sleptons, e.g., X2 — * t + t, f — > z/ + / + /'. 

We see from Fig.[2]that there are several possible scenarios for the sparticle mass spectrum 
in the NUHM, which are distinct from the standard CMSSM spectrum as we now enumerate. 

1. > rngj^ > rrii^^ > ruf > rriy^: 

This is the mass hierarchy seen in panel (d) of Fig. [2] for larger itla- The neutralino 
can decay into 

Z>e + Ve 

Ti + r 

Vr + Vt (15) 
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which, for the first mode, would be followed by 



the second mode by 



~eL ^ i>e + f + / 

fi + e + r 
f 1 + z/e + z/^ 

h'r + l^e + T, (16) 



i>e fi + Ue + T 

fi + e + Ur 

Vr + Ve + 

i>r + e + T, (17) 



and the third mode by 

n^i>r + f + f 

In the case of the two-body neutralino decay to stau and tau, the decay rate is 



. , + + mi — 2{m'im% + mlml + mlml) 

ix-^n + T) = ^ — — 

32TTmi 

A 

X {i\Cn\' + \CL\')iml + ml- ml) - 2(Cz.Q + CnCl)m,m, 

(19) 

where Cl, Cr are the left and right couplings in the neutralino-tau-stau vertex. There 
are similar expressions are for the other two-body decay modes. 

2. m^ > mo^ > mf > mg^ > mc,^: 

This hierarchy occurs for more intermediate values of m^ when is large as seen in 
panels (b) and (c) of Fig. The neutralino 2-body decay modes would be the same 
as in the previous case, although with different branching ratios. However, the cascade 
decays are in general different. In this case, we would have 

^ h + f' + f 



^•^The viability of such models would require some action to conform with the GUT constraint, e.g., by 
constraining inflationary cosmology. 
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eL + e + Ur 

ei + Ve + T 
i)e + Ue + 
Ur + e + T 



(20) 



f eL + e + T 

eL + l^e + l^T 
Ue + l^e + T 

D^ + e + Vr (21) 

and 

~eL^i^e + f + /. (22) 

3. For large tan/3, e.g. tan/3 = 40 as shown in panel d of Fig. [3], and large the 
third-generation sleptons get larger masses through the Yukawa couplings. Thus we 
have > m^^ > m^^ for the lightest sparticles. In this case, the neutralino cascade 
decays become simpler, 

Ve + (23) 

which, for the first mode, would be followed by 

iL^-Ve^ f + /, (24) 

whilst the second mode is invisible. 

4. There are also other possibilities for narrower region of parameter space, near where 
the masses cross each other in Fig. [2J (a) m,>^ > rrif > rngj^ > > nii^^; (b) 

> nif > > nie^ > rriy^] and (c) rriy,^ > > rrif > rrigj^ > rriy^. 

We extract from these examples a few generic features. As in the cases of many other 
scenarios beyond the Standard Model, particularly within the general framework of super- 
symmetry, the most prominent signature of a sneutrino NLSP scenario is likely to be missing 
energy. However, there would in general be accompanying signatures that would enable a 
sneutrino NLSP scenario to be distinguished from other possibilities. Specifically, one ex- 
pects to see also events with missing energy accompanied by leptons. The precise nature of 
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this supplementary signature would, however, depend on the nature of the sneutrino: z/^, 
or i>e, and on the hierarchy of heavier sparticle masses. 

In particular, the relative multiplicities of different charged leptons would depend on the 
flavour of the invisible and hence be a useful tool for identifying it. Concretely, in cases 
where the parent sparticle has no lepton flavour, as would normally be the case at the LHC, 
each Vt NLSP would be accompanied by an unmatched r or i/^, and each z/e^^u NLSP would 
be accompanied by an unmatched electron, // or corresponding neutrino. In general, there 
would be additional lepton- antilepton pairs with matched flavours. 

7 Summary 

We have analyzed in this paper the possibihty of a sneutrino NLSP in NUHM models with 
a gravitino LSP. This possibihty does not exist in the CMSSM, but is quite generic in the 
NUHM, as we have illustrated with various speciflc examples. The sneutrino might well be 
the Vr, but the and z/g are also possible candidates for the NLSP. A sneutrino NLSP would 
be metastable and subject to cosmological constraints on late-decaying particles, but we have 
shown that these are not difficult to respect. There are various different possible scenarios for 
the spectrum of sparticles heavier than the sneutrino, which would have distinctive signatures 
at colliders. In addition to events with missing energy carried away by the invisible z?, there 
would also be events with accompanying charged leptons. The flavours of such leptons would 
help identify the flavour of the sneutrino NLSP. 

As particle physics embarks on the study of the TeV scale with the LHC, much unknown 
physics will surely be revealed. Supersymmetry is occasionally regarded as a 'known un- 
known' in the sense that, whereas we do not know whether it exists,we think we know what 
it would look like if it does exist. This paper reminds us that supersymmetry should rather 
be regarded as an 'unknown unknown', in the sense that not only do we not know whether 
it exists, but we also do not know what it would look like. In the conventional 'known 
unknown' scenario, the LSP is the lightest neutralino and supersymmetry would produce 
missing-energy events. The latter would also be the signature of a scenario with a gravitino 
LSP with a neutralino NLSP, at least in gravity-mediated scenarios. However, once the 
'Pandora's box' of a gravitino LSP has been opened, many other NLSP candidates fly out. 
In addition to the relatively familiar case of the lighter stau and the more radical case of the 
lighter stop, there are other possibihties including the sneutrino NLSP scenarios discussed 
here. All of these scenarios have distinctive features, as illustrated here, so the LHC and 
subsequent experiments have good prospects for detecting and distinguishing between the 
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various 'unknown unknowns'. No 'unknown unknown' stone should be left unturned in the 
search for supersymmetry. 
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Appendix: Sneutrino Three-Body Decay 

We calculate here the radiative sneutrino three-body decay 



(25) 



that may arise through the photino content of the neutralino, as illustrated in the diagram 
below. 




The invariant amplitude for this decay is 



M 



4Mpi(g2 - mj 



(26) 



where M-p\ = 1I\/WkGn is the Planck mass, rrii = m^o, 



q = k + Pi 



G 



and the dimensionless couplings are 



B, 



92 

V2 



{02i - tanOwOii) 



d = Oii cos ^vi/ + sin 



(27) 

(28) 
(29) 
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Ignoring the neutrino mass, we get 



-72 Sk- Pg 



^3 ^i^j 



G 



2(Pg ■q){p-q) + {Pg ■ p){mimj - q ) 



+2{k ■ q){p- q) + {k ■ p){mimj - g^)} ■ 



(30) 



Use the Dahtz parametrization mu = pi+ p2, = P2+P3 with pi — Pq, P2 — k, ps — p, 
we get 



Pg 


k = 




(31) 


k 


P = 




(32) 


Pg 


P = 




(33) 


k 


Q = 




(34) 


P 


Q = 




(35) 


Pg 


Q = 


|(m?2 + "^g) 


(36) 



and 



2 2 
q =mi2. 



(37) 



Hereafter, we abbreviate our notation by defining M = mo and ttig = mg. 
The resulting partial decay rate is [37] 



1 



(38) 



(27r)3 32M3 ' 

This can be integrated analytically using the following integration boundaries - for ^33: 
and (771^2^2 + m^2"^G ~ "^12 ~ -^^^g)/"^i2) for ml2- ^Uq to M^. The result is 



7687r3M|iM3 



(39) 



where 



ml + m| 
mQ{mimjY 
1 

mQimimjY 



X:^^(M^«-m^") + ao(i,j)ln 



.a=l 



M2 



^«a+iM(M2«-mg') + ai(i,j)ln 

.a=l 

1 ^ 



G 

M2 



(40) 



Gj 



-Q;o(«,i) 



(m2 
23 



(41) 



+/?o(^,i)ln 



G 



2 



and 



Id{i,j) = Ic{j,i), 

where the auxihary functions a and P are defined below. Note that there is 
singularity when i — j, because Ic + Id is of the form 

ff{a,b) f{b,a)\ 1 
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a — b 



(42) 

(43) 
no actual 

(44) 



and in this case, with rrii = mj, we get 
Ic + Id ^ ^ 



2 2 



-2— + (-13m; - 2^6) — + (-35m| - Umfae - Sag) — 
I u o 



u 



+(-49m^ - 24m^Q;6 - 9m- as - 2q;4) — 



+ (— 35m^ — 25m^a6 — ISm^as — 7mj^a4 — 203) — 



+(— 7mj — lOmjCte — lOm^as — 8m^a4 — 5m^as — 2a2) — 
+ (7mJ^ + 3m]^°Q;6 — 2m^Q;4 — 3m^Q;3 — "iw^^a^ — 2ai)M 
+ (5mJ^ + Aml^aQ + Sm^^^as + 2m^a4 + m^as — m^ai — 2ao) Inu 



.10 



1 



u. 



(45) 



u=rrfp^—m^ 



The auxiliary functions are 

ai{ij) = mimj{8m%M'' + 6M'^m%) - 3m%M^ 
"2(^,7) 
"3(^,7) 
Q;4(i,j) 
a5(«, j) 
a6(«,i) 

a7(^,i) = 1, 



G 



mirrij [ om, 



3m 



G 



mimj{M* - 6m%) + 8m% + 12M^m% 

- 6ma - 2mimjM'^ 
mirrij - 2M^ 



2_4 



(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
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and 



(54) 



a=0 



E2(a-1) /. .X 



(55) 



a=l 



/52(^,j) 

PS. 3) 

PeiiJ) 



35m^ + 20m^aQ + lOm^a^ + Am^a^ + 
35mf + 15mjaQ + bm'^a^ + 

7m^ + aQ 



(56) 
(57) 
(58) 
(59) 
(60) 
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